Ferromagnetism in a lattice of Bose condensates 
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We show that an ensemble of spinor Bose-Einstein condensates confined in a one dimensional 
optical lattice can undergo a ferromagnetic phase transition and spontaneous magnetization arises 
due to the magnetic dipole-dipole interaction. This phenomenon is analogous to ferromagnetism in 
solid state physics, but occurs with bosons instead of fermions. 
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The Heisenberg model of spin-spin interactions is con- 
sidered as the starting point for understanding many 
complex magnetic structures in solids. In particular, it 
explains the existence of ferromagnetism and antiferro- 
magnetism at temperatures below the Curie tempera- 
ture. It is defined by the spin Hamiltonian Q 



H, 



spin ^ Jij S j • Sj , 



where Si is the spin operator for i-th electron, and are 
known as the exchange coupling constants. This Hamil- 
tonian arises from the direct Coulomb interaction among 
electrons and the Pauli exclusion principle. In addition to 
the exchange interaction, there exists another important 
type of magnetic interaction, the magnetic dipole-dipole 
interaction. However, in solid materials the dipolar cou- 
pling is typically several orders of magnitude weaker than 
the exchange coupling, and would correspond to Curie 
temperatures much below the observed ones. Hence its 
contribution to the spin Hamiltonian can be neglected 
in practice. Q It follows fr om this argument that ferro- 
magnetism is not generally expected to occur in bosonic 
lattices of neutral atoms, a result of the inapplicability of 
the Pauli principle, the absence of Coulomb interaction 
and small atomic magnetic dipole moments. However, 
we qualify this remark by noting that as a result of ac- 
curate recent measurements of the scattering length of 
spin-changing collisions, it is now established that the 
ground state of optically trapped 87 Rb spinor Bose con- 
densates is ferromagnetic, an important result that we 
use in the following |^,||. 

The goal of this paper is to show that this result, com- 
bined with the recent experimental realization of regular 
arrays of Bose-Einstein condensates in optical lattices, 
leads to a situation where it becomes possible to carry 
out detailed static and dynamic studies of magnetism on 
one to three-dimensional periodic lattices. 

Quantum degenerate Bose gases on optical lattices 
were first used to demonstrate "mode-locked" atom lasers 



1 We note however that the magnetic dipole-dipole interac- 
tion plays an important role in domain formation in macro- 
scopic samples. 



H . It has also been theoretically demonstrated that they 
undergo a Mott insulator phase transition as the depth 
of the lattice wells is increased || . Recently, they have 
become of interest in the study of quantum chaos 
Here, we show that spinor condensates, localized in opti- 
cal lattices deep enough for the individual sites to be in- 
dependent, can undergo a ferromagnetic-like phase tran- 
sition that leads to a "macroscopic" magnetization of the 
condensate array. 

We consider specifically the case of spinor 87 Rb con- 
densates which are as we have mentioned individual 
ferromagnets of random directions in the absence of ex- 
ternal fields and magnetic dipole-dipole interaction. We 
show that the magnetic dipole-dipole interaction between 
lattice sites can spontaneously align the magnetization 
of the individual sites. This is possible because of the 
Bose enhanced magnetic dipole moments of the conden- 
sate which in turn enhances the strength of the magnetic 
dipolar interaction. 

Our starting point is the Hamiltonian H describing an 
F = 1 spinor condensate at zero temperature trapped 
in an optical lattice, subject to a magnetic dipole-dipole 
interaction H^d and is coupled to an external magnetic 
field via the magnetic dipole Hamiltonian Hb PTES], 



H = H Q + H t 



dd 



(1) 



Here 



H = J rf 3 rV4(r)h?i 2 V 2 /2m + ^(r)]^(r) 
+ (A s /2) j d 3 r^ a (r)ijl(r)vj (r)yj a (r) 
+ (A Q /2) / d 3 r^(r)^(r)F Q/3 .F^(r)^(r) 



describes the interaction of the atoms with the lattice 
potential Vj,(r) and ground state collisions. It includes 
an implicit sum over the indices {a, (3, fj,, v} = {0, ±1} 
that label the three Zeeman sublevels. The parameters 
A s and A a characterize the short-range spin-independent 
and spin-changing s-wave collisions, respectively. Specif- 
ically, A a is proportional to the difference between the 
s-wave scattering lengths in the triplet and singlet chan- 
nels jl0| |l2| . It has recently been measured to be negative 



1 



for Rb, accounting for its ferromagnetic ground state where \' a 



Our model includes the long-range magnetic dipole- 
dipole interaction between different lattice sites, but ne- 
glect it within each site, assuming that it is much weaker 
than the s-wave interaction described by Hq. We also 
assume that the optical lattice potential is deep enough 
that there is no spatial overlap between the condensates 
at different lattice sites. We can then expand the atomic 
field operator as ip(r) = HiEa^o ±1 *a(')4( r ) where i 
labels the lattice sites. 

The Hartree wave function 4>i(v), determined by min- 
imizing the total energy, is the wave function of the con- 
densate at the i-th site. It is assumed that all Zeeman 
components share the same spatial wave function. If the 
condensates at each lattice sites contain the same num- 
ber N of atoms, then the ground-state wave functions 
for different sites have the same form (f>i(r) = cf>(r — rj). 
Under this condition, the dipolar interaction potential is 
Had 



= EijYi V dd with 



y dd — 



Mo 
4tt 



(2) 



where po is the vacuum permeability, = — ij, 
r.jj = Tij/\rij\, r-j is the coordinate of the i-th site, and 
Pi = "/Si is the magnetic dipole moment at site i, with 
Sj = a\ c (i)~F a pap{i) being the total angular momentum 
operator and 7 = gFPB the gyromagnetic ratio. Finally, 
the coupling of the atoms to the external magnetic field 
B cxt is described by 



H, 



= - 7 ^S, -B c 



In this letter we consider a one-dimensional optical lat- 
tice along the z-direction, which we also choose as the 
quantization axis. Hence 



V. 



y _ 

dd 



4 7 r|r lJ | 3 



[fli ■ ftj - Zp z t pf\ , 



with pA being the z-component of pj. In the absence 
of long-range magnetic dipole-dipole interaction and of 
external magnetic fields, the individual condensates can 
therefore be considered as independent "magnets" whose 
pseudo-spin vectors point in random directions, with no 
spin correlations between sites. Our goal is to determine 
the spin structure of the system if the different sites are 
allowed to interact with each other through the magnetic 
dipole-dipole interaction. 

In the absence of spatial overlap between individual 
condensates, and neglecting unimportant constants, the 
total Hamiltonian of the system takes the form p2| 



(3) 



(l/2)A a /d 3 r|^(r)| 4 and A, 



7Mo/(47r|r y f). 

In general the external magnetic field consists of two 
contributions: a controlled, external applied field; and 
an effective "stray" field that accounts for all possible ef- 
fects from the experimental environment and the system 
itself. A typical example is the environmental magnetic 
fluctuations. In the present case, we take the applied field 
along the quantization axis z. The effective environmen- 
tal magnetic field, in contrast, can have any orientation. 
We choose it without loss of generality to be along the 
transverse direction, including any longitudinal compo- 
nent in the definition of the applied field. Hence the 
external field has the form 



B, 



B z z + B p p, 



where p = y/x 2 + y 2 is the radial coordinate. 

If the optical lattice is sufficiently long, one can safely 
neglect the boundary effects and concentrate on a generic 
single site i of spin S. Its Hamiltonian reads 



h = \' a S 2 - 1 S 



B z 



2$> 



B n 



(4) 



We determine its ground state in the mean-field ap- 
proximation (also known as the Weiss molecular poten- 
tial approximation) |l3|] . It consists in replacing the op- 
erators S" , a = p, z, by their ground-state expectation 
value 



(Sf) 



Nm n 



(5) 



which is assumed to be the same for different sites. We re- 
mark that m z is nothing but the difference in population 
of the Zeeman sublevels of magnetic quantum numbers 
±1. This allows us to approximate the Hamiltonian (Q) 
by 



h mi = \' a S 2 — 7S • B 



eff , 



(6) 



where we have introduced the effective magnetic field 
B c ff = (B z + 2Km z )i + (B x - Am p )p 



(7) 



andA = iV£ j¥i Aij. 

We mentioned that the individual spinor condensates 
at the lattice sites are ferromagnetic, < 0. In that 
case, the ground state of the mean-field Hamiltonian (||) 
must correspond to a situation where the condensate at 
the site i under consideration must be aligned along B e g 
and takes its maximum possible value N. That is, the 
ground state of the mean-field Hamiltonian (^|) is simply 



\GS) = \N,N) BaH , 



(8) 
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where the first number denotes the total angular momen- 
tum and the second its component along the direction of 
Beg. Note that \GS) represents a spin coherent state 
in the basis of \S, S z ). The fact that the ground state 
magnetic dipole moment of each lattice site is N times 
that of an individual atom results in a significant mag- 
netic dipole-dipole interaction even for lattice points sep- 
arated by hundreds of nanometers. This feature, which 
can be interpreted as a signature of Bose enhancement, 
is in stark contrast with usual ferromagnetism, where the 
magnetic interaction is negligible compared to exchange 
and where the use of fermions is essential [[l3| . 

The mean- field ground state of Eq. (|^) allows us to cal- 
culate the magnetization components m z and m x . One 
finds readily 



1 

N 



(GS\S?\GS) =cos( 



(9) 



where 8 a is the angle between B c g and the a-axis. 
For B z = 0, Eq. (|) yields 

2Am, 



^(2Am 2 ) 2 + (B p - Am p f 
B p - Am p 



(10a) 
(10b) 



9 y/{2h.m z y + {B p -hm p y 
The solutions to Eqs. ( |Io| ) can be divided into two cases: 

1. For B p > 3 A, the only solutions are m z — and 
772 p — 1. That is, the lattice of condensates is mag- 
netically polarized along the transverse magnetic 
field. 

2. For Bp < 3 A, there are two coexisting sets of 



solutions: (i) 777 z = and m p 



1 ; and (ii) 



777 z = ±yjl - (B p /3A) 2 and m p = B p /3A. How- 
ever, it is easily seen that the state associated with 
the latter solutions has the lower energy. Hence it 
corresponds to the true ground state, while solution 
i) represents an unstable equilibrium. 

We have, then, the following situation: As the environ- 
mental effective magnetic field strength is reduced below 
a critical value 3A, the lattice ceases to be polarized along 
the direction of that field. A phase transition occurs, 
and a spontaneous magnetization along the z-direction 
appears, characterized by a finite m z . 

This phenomenon is reminiscent of conventional ferro- 
magnetism. Indeed, our model is somewhat analogous 
to the Ising model Q , with the environmental transverse 
magnetic field B p playing the role of temperature. For 
Bp = — corresponding to zero temperature in Ising 
model — the spins at each lattice site Si align them- 
selves along the lattice direction, even in the absence of 
longitudinal field. This spontaneous spin magnetization 
diminishes as B p increases, and completely vanishes if B p 
exceeds the critical value 3A — the analog of the Curie 
temperature in the Ising model. We note however that 



the two cases exhibit important qualitative differences: 
For example, no spontaneous magnetization occurs in ID 
Ising model, for any finite temperature |Q . 

Our analysis so far assumes an infinite one-dimensional 
lattice. This assumption is required for the validity of 
mean-field approximation. We note however that the ap- 
pearance of a spontaneous magnetization does not rely on 
this condition being fulfilled. We demonstrate this point 
by considering just two lattice sites, that is, a double- 
well situation. In this case, we can numerically solve 
the Hamiltonian (^) , without invoking the mean- field ap- 
proximation, by expanding the Hamiltonian matrix on 
the basis oi\S x = N, Sf ) ® \S 2 = N, S| ). Figure |l| shows 
the spontaneous magnetization m z as a function of the 
external field strength B p . As the Hamiltonian matrix 
has a size of (2A^ + l) 2 x (2N + l) 2 , we cannot use a very 
large N. However, as we can see from the figure, as N 
increases, m z rapidly approaches the mean-field results. 
Hence the mean-field treatment is in fact a surprisingly 
good approximation for our purpose. 




Bp /A 

FIG. 1. Spontaneous magnetization as a function of envi- 
ronmental magnetic field strength. Only the positive values 
are plott ed. The dashe d line represents the mean-field result 
m 2 = yjl — (B X /3A) 2 and the solid lines, from left to right, 
correspond to the exact numerical results for a two-site lattice 
with N = 10, 15 and 25 atoms, respectively. 



To estimate the feasibility of an experimental detec- 
tion of the spontaneous magnetization, we consider as 
an example the F = 1 electronic ground state ( 3 Si/2) 
of 87 Rb, for which A' < The Lande factor for 

this state is = —1/2. For an optical potential of pe- 
riod equal to 426 nm (nearest neighbor separation) we 

find A = NJ2i^i X v ~ 1 - m x 10 " 7 G - If the Particle 
number at each site is N = 2000, this gives a critical 
value of the environmental magnetic field of 3A = lmG. 
These numbers indicate that the observation of sponta- 
neous magnetization in a lattice of spinor condensates is 
well within experimental reach. 

In a typical experimental situation, the environmental 
magnetic field is likely to exhibit temporal fluctuations. 
To account for them, we assume that this field has a 
uniform angular distribution and that the fluctuations in 
field strength have a width A^. The time-averaged en- 
vironmental field is again assumed to be along the trans- 
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verse direction, with strength B p . Assuming the magne- 
tization of the system follows the field fluctuations adia- 
batically, it is not difficult to include this effect into our 
mean-field treatment. The time-averaged spontaneous 
magnetization along z-axis under satisfies then the cubic 
equation 

(1 - m2 )(9A 2 m 2 + A | /2) = {B 2 + A 2 B)ml 

Figure || illustrates m z as a function of B p for various 
As- As we can see, the effect of the field fluctuation 
is to decrease the maximum spontaneous magnetization, 
and to increase the critical field strength at which the 
spontaneous magnetization vanishes. The maximun m z 
at B p — decreases from 1 to 1/3 when Ab increases 
from to infinity. 



is anisotropic and contains both ferromagnetic and anti- 
ferromagnetic terms. In higher-dimensional lattices, the 
ground-state spin structure will thus become much richer. 

We conclude by remarking that in addition to this in- 
trinsic interest, the study of the magnetic properties of 
spinor condensate lattices provide us with a highly con- 
trollable test system to study fundamental static and 
dynamical aspects of magnetism and lattice dynamics, 
including, e.g., the role of dimensionality in phase tran- 
sitions and macroscopic quantum tunneling of the mag- 
netic moments. 

This work is supported in part by the US Office of 
Naval Research under Contract No. 14-91-J1205, by the 
National Science Foundation under Grant No. PHY98- 
01099, by the US Army Research Office, by NASA, and 
by the Joint Services Optics Program. 
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Bp /A 

FIG. 2. Mean-field results of the time-averaged sponta- 
neous magnetization as a function of environmental magnetic 
field strength, taking the magnetic field fluctuations into ac- 
cout. Only the positive values are plotted. 



Recently, condensates trapped in periodic optical po- 
tentials have attracted much attention. The phase coher- 
ence and atom statistics of the system have been studied 
experimentally |^,|l^,^5). However their magnetic prop- 
erties have not been fully explored. A few years ago, 
Meacher et al. observed experimentally the paramag- 
netic behavior of cold (but non-condensed) cesium atoms 
confined in an optical lattice jl6| . Here we have demon- 
strated that by replacing the cold atoms with a spinor 
condensate, ferromagnetism can be observed. This is 
made possible by the collectively enhanced magnetic mo- 
ments of the condensate which in turn enhances magnetic 
dipole-dipole interaction between different lattice sites. 

In future work it will be interesting to study the dy- 
namical response of the system under the effect of an ex- 
ternal time-dependent longitudinal magnetic field. The 
instantaneous magnetization in this case is likely to form 
hysteresis loops which might find applications e.g. in 
quantum information processing. Further studies should 
also include the properties of finite temperature exci- 
tations of the system, which correspond to spin waves. 
The extension of this work beyond one-dimensional lat- 
tices also shows much promise: Although it is predomi- 
nantly ferromagnetic in ID, the dipole-dipole interaction 
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